Dynamical slowdown of polymers in laminar and random flows 
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The influence of an external flow on the relaxation dynamics of a single polymer is investigated 
theoretically and numerically. We show that a pronounced dynamical slowdown occurs in the vicinity 
of the coil-stretch transition, especially when the dependence on polymer conformation of the drag 
is accounted for. For the elongational flow, relaxation times are exceedingly larger than the Zimm 
relaxation time, resulting in the observation of conformation hysteresis. For random smooth flows 
hysteresis is not present. Yet, relaxation dynamics is signiflcantly slowed down because of the large 
variety of accessible polymer configurations. The implications of these results for the modeling of 
dilute polymer solutions in turbulent flows are addressed. 
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The dynamics of an isolated polymer in a flow field 
forms the basis of constitutive models for dilute polymer 
solutions 0,0,1^' The modeling of drag- reducing flows, 
for instance, requires an appropriate description of sin- 
gle polymer deformation in turbulent velocity fields 0. 
In the last decade, major advances in fiuorescence mi- 
croscopy offered the possibility of tracking isolated poly- 
mers both in laminar Q and random flows 0. The dy- 
namics of a polymer in thermal equilibrium with the sur- 
rounding solvent is commonly described in terms of nor- 
mal modes and relaxation times associated with them. 
The analytical form of the relaxation spectrum was first 
obtained by Rouse under the assumption that the poly- 
mer could be described as a series of beads connected by 
Hookean springs . The Rouse model was subsequently 
improved by Zimm to include hydrodynamical interac- 
tions between segments of the polymer 7:]. In Zimm's 
formulation, the equations of motion are decoupled into 
a normal mode structure by preaveraging the distances 
between the beads over the distribution of polymer con- 
figurations. The relaxation time associated with the fun- 
damental mode, T, determines the typical time that it 
takes for a deformed polymer to recover the equilibrium 
configuration in a solvent. The normal mode theory was 
confirmed by the analysis of the oscillatory motion of a 
DNA molecule immersed in a solvent and held in a par- 
tially extended state by means of optical tweezers 8j . An 
alternative approach to examine polymer relaxation con- 
sists in stretching a tethered DNA molecule in a fiow and 
measuring its relaxation after cessation of the fiow [gllioj. 
The theoretical predictions for these experimental con- 
ditions are provided by the scaling theory and the 
static dynamics formalism jl^ . 

The aforementioned studies all consider the internal 
dynamics of a polymer floating in a solvent under the 
influence of thermal noise only — the interaction of the 
polymer with an external flow is not taken into account. 
One of the aspects highlighted by experiments is that 



polymer dynamics in a moving fluid is strongly influenced 
by the carrier velocity field. The coil-stretch transition is 
the most noticeable example: as the strain rate exceeds a 
threshold value, the polymer undergoes a transition from 
the coiled, equilibrium configuration to an almost fully 
extended one [ll^ . Therefore, when a polymer is freely 
transported by a non-uniform fiow, we expect that the 
time scales describing its dynamics may be significantly 
different from the Zimm time r. Simple models for poly- 
mer stretching indeed suggest deviations from Zimm's 
theory near the coil-stretch transition |illil[l3. Dis- 
crepancies in the definition of the correct relaxation time 
are also encountered in drop formation experiments |l7j . 

In this Letter we investigate polymer relaxation dy- 
namics both in elongational and random smooth flows. 
Our analysis brings to evidence an important slowdown 
of dynamics with respect to the Zimm timescales, in 
the vicinity of the coil-stretch transition. For the elon- 
gational fiow, this is related to conformation hystere- 
sis 0, 0, 113. For random fiows, we show that hys- 
teresis is not present. Nonetheless, the amplification of 
the relaxation time persists, albeit to a lesser extent, due 
to the large heterogeneity of polymer configurations. In 
both cases, the dependence of the drag force on the poly- 
mer configuration plays a prominent role. This suggests 
the necessity of improving current models of polymer so- 
lutions in turbulent fiows to account for such effect. 

The dumbbell approximation is the basis of the most 
common models of single polymer dynamics and vis- 
coelastic models of dilute polymer solutions jl^l. Its va- 
lidity relies on the fact that the slowest deformational 
mode of the polymer is the most influential in produc- 
ing viscoelasticity 0. However, when attention is di- 
rected to non-equilibrium dynamics it is often too crude 
to assume that r is independent of the conformation of 
the molecule [i^. Theref ore, following de Gennes 0|, 
Hinch [2^, and Tanner's |2J| approach, we consider a 
model where the polymer is described as two beads con- 
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nected by an elastic spring and the probability density 
function of the end-to-end vector, ^'(i?, t), satisfies the 
diffusion equation: 
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where Hijit) — djVi{t) is the velocity gradient, Rq is the 
mean extension at thermal equilibrium and R — \R\. 
The function f{R) defines the entropic force restoring 
stretched molecules into the coiled configuration. Syn- 
thetic polymers are properly described by the Warner 
law, f{R) = 1/(1 — -R^/L^), where L is the contour 
length of the polymer; biological macromolecules are 
better characterized by the Marko-Siggia law, f{R) = 
2/3-L/{6R) + L/[6R{l-R/L)^] The fiow strength 
relative to the polymer tendency to recoil is expressed 
by the Weissenberg number Wi, defined as the product 
of the Zimm time r by a characteristic extension rate 
of the flow. The function encodes for the depen- 

dence on the polymer conformation of the drag exerted 
by the fluid: a spherical coil offers a smaller resistance 
with respect to a long rod-like configuration. We utilize 
the expression i'{R) = 1 -|- (Cs/Cc — 1)R/L that interpo- 
lates linearly between these extremes (see Refs. (25tl2q). 
Here, Cc = SVOt^Roti/S and (s — 271 Lt] / ln{L / i) are the 
friction coefficients for the coiled and the stretched con- 
figuration, respectively, 77 is the solvent viscosity and £ 
is the diameter of the molecule. Strictly speaking, the 
above form of viR) was deduced from experiments and 
microscopic simulations of laminar fiows. To our knowl- 
edge, measures of J^(i?) in random flows are not available 
yet, and the study of its functional dependence lies be- 
yond the scope of the present work. Assuming that the 
polymer is aligned with the local stretching direction of 
the flow for the most part of its evolution, we shall use a 
linear ^{R) for a random flow as well. 

To define the relaxation time in the presence of an ar- 
bitrary external fiow, we consider the probability density 
function of the rescaled extension, P(r, t) with r = R/L, 
and identify the dynamical relaxation time, U-d, as the 
characteristic time needed for P(r, t) to attain its station- 
ary form Pstir). In the cases examined in this Letter, as 
we shall see, the probability density function of r satisfies 
the Fokker-Planck equation 
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where the form of -Di(r) and D2{r) depends on the 
flow and t' — t/r. The stationary solution to Eq. (j^J 
takes the potential form Pst{r) = N e-yiY>[—E{r) / KbT], 
where N is the normalization constant and E{r) = 
—KbT Di{p) / D2{p) dp. The finite-time solution ad- 
mits the expansion 
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P{r, t') = Pst (r) + a„p„(r) e"*'/-- , 



(3) 



Figure 1: Elongational flow: (a) rescaled relaxation time 
vs. Wi for b = 400. The entropic force is given by the Warner 
law. For small Wi, ty-Bi/r grows as 1/(1 — 2 Wi), as can be seen 
by replacing /(r) with 1 in Eq. Q- For Wi » Wicrit, the time 
needed to reach the asymptotic regime is set by the time scale 
of the flow, 7~^, and trei/T decreases as Wi~^ . (b) Rescaled 
maximum relaxation time tmax/T vs. Cs/Cc (6 = 400). 



where the coefficients a„ are fixed by P(r, 0), Pn{f) are 
the eigenfunctions of the Fokker-Planck operator, and (t„ 
are the reciprocals of its (strictly positive) eigenvalues, 
arranged in descending order (cr„ > Un+i)- The relax- 
ation time is thus deflned as iici = uit. 

As a first example, we examine the steady planar elon- 
gational flow V = 7(0;, —y). By assuming that the poly- 
mer extension in the x-direction is much greater than in 
the ?/-direction, it is easy to derive from Eq. (^) an equa- 
tion of the form Q with P>i(r) ^ Wir - f{r)r/[2v{r)], 
D2{r) = [2bv{r)]-\ f{r) = f{rL), v{r) = v{rL), 
b = L^/Rl, and Wi = 7T For (s = Cc ^rci can be 

computed from by solving a central two-point con- 
nection problem for a generalized spheroidal wave equa- 
tion In the general case, Cs > Cc, we resorted to a 
numerical computation based on the variation-iteration 
method of quantum mechanics l28j. In the vicinity of 
the coil-stretch transition {Wi — 1/2) shows a sharp 
peak (Fig. In this range of Wi there is a critical 
competition between the entropic force and the velocity 
gradient that makes the convergence time to the steady 
state extremely long. This effect is strongly enhanced 
by the conformation-dependent drag; the peak tmax in- 
deed increases with Cs/Cc (Fig.QJ. Those extremely long 
relaxation times arc intimately connected with the obser- 
vation of flnitc-time conformation hysteresis ^18, .19. 2(|. 
For large enough Cs/Cc there is a narrow range of Wi 
around the critical value for the coil-stretch transition 
where E{r) has a double well structure ilSt^, JCy. The 
barrier height separating the coiled and the stretched 
state is much greater than the thermal energy KbT, and 
therefore the polymer remains trapped in its initial con- 
flguration for an exceptionally long time (Fig. ^ . 

The discovery of elastic turbulence has recently al- 
lowed the examination of single polymer dynamics in a 
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Figure 2: (a) Effective energy at the coil-stretcli transition for 
a polyacrylamide (PAM) molecule {b = 3953, Cs/Cc = 6.87) 
in the elongational flow (dashed line) and the 3D Batchelor- 
Kraichnan flow (solid line). The left vertical axis refers to 
the dashed line; the right vertical axis refers to the solid line; 
(b) effective energy in the random flow for a PAM molecule 
with constant drag (Cs=Cc) (from top to bottom Wi — 
0.4,0.7,1.0,1.3); (c) the same as (b), but with = 6.87Cc 
(from top to bottom = 0.28, 0.33, 0.38, 0.43). When > 
C^c the transition occurs in a much narrower range of Wi. 



random smooth flow generated by viscoelastic instabili- 
ties 0, . The velocity gradient fluctuates along fluid 
trajectories; the average local deformation rates deflne 
the three Lyapunov exponents of the flow. Experimental 
observations have been accompanied by theoretical and 
numerical studies based on the dumbbell model (see, 
e.g., Refs. 0, 0, H^l)- To analytically investigate 
polymer relaxation dynamics in random flows, we ini- 
tially assume that the velocity field obeys the Batchelor- 
Kraichnan statistics [sJl . The velocity gradient is then a 
statistically isotropic and parity invariant Gaussian pro- 
cess with zero mean and correlation: {K,ij{t)nki{s)) = 
2\5{t-s)[{d+\)5ikSji-5ij5ki-5ii5jk\l[d{d-l)], where dis 
the dimension of the flow and A denotes the largest Lya- 
punov exponent. In this context, we indicate by P(r, t') 
the probability density function of the extension both 
with respect to thermal noise and the realizations of the 
velocity fleld. For the elongational flow Eq. (|2Jl was ob- 
tained under the uniaxial approximation. On the con- 
trary, for the isotropic random flow Eq. |(5J holds ex- 
actly and can be obtained by a Gaussian integration by 
parts followed by integration over angular variables. The 
drift and diffusion coefficients take the form Di{r) = 
{d - l)/d Wi r - f{r)rl[2v{r)] + {d - l)/[2bi){r)r] and 
D2{r) = Wt r^/d+ [25i>(r)]"i with Wi = At. The sta- 
tionary distribution admits once more a potential form. 
For d = 3 and Cs > Cc the potential displays a very wide 
well, the effect of the conformation-dependent drag being 
to increase the probability of large extensions. There is 
no evidence of pronounced double wells (Fig. EJ. Near 
the coil-stretch transition, the effective barrier heights 
separating the coiled and stretched states are indeed at 



most comparable to thermal energy. For realistic Cs/Cc 
no conformation hysteresis is therefore expected to be 
observed in random flows. The behavior of tj-ci vs. Wi 
is however analogous to the one encountered in the elon- 
gational flow: tj-ei/r increases as [1 — Wi{d + 2)/d]~^ at 
small Wi and decreases as Wi~^ at large Wi. A peak 
near the transition is present, that becomes more and 
more pronounced with increasing Cs/Cc, attaining values 
as large as about thirty times r (Fig. . The reason for 
this behavior is the breadth of Pst {r) and the consequent 
large heterogeneity of accessible polymer configurations. 

To corroborate the results obtained in the context of 
the short-correlated flow, we performed Brownian Dy- 
namics simulations of dumbbell molecules [s^ in the ran- 
dom flow introduced by Brunk et al. |34l |. This model re- 
produces the small-scale structure of a turbulent flow by 
means of a statistically isotropic Gaussian velocity gra- 
dient. The autocorrelation times of components of the 
strain and rotation tensors are set to be multiple of the 
Kolmogorov time r,, by comparison with direct numer- 
ical simulations of 3D isotropic turbulence {ts = 2.3r^, 
tr = 7.2t^). The Lyapunov exponent of this flow is 
A ~ lOr^^. We computed irei as the time of convergence 
of the moments (r"(t)) to their stationary value (r")st: 
t^el = -hmt^ooln[(r"(t)) - (r")st]/t, where the aver- 
ages were taken over an ensemble of realizations of the 
flow, all with the same initial extension r(0). The nu- 
merical difficulty arising from the singularity of the en- 
tropic force at R = L has been overcome by exploiting 
the algorithm introduced in |35j . The results shown in 
Fig.Elconfirm the scenario depicted in the context of the 
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Figure 3: 3D Batchelor-Kraichnan flow: (a) trd/r vs. Wi 
for a PAM molecule (h = 3953); (b) tmax/T for the following 
polymers: DNA (•, h = 191.5; o,h = 260; □, fo = 565; +,h^ 
2250), polystyrene (x, fe = 673), polyethyleneoxide (PEO) 
(A, & = 1666), Escherichia Coh DNA (A, & = 9250), PAM 
(■). Measures of b and Cs/Cc can be found in 18, 19, 20, 25J. 
Synthetic polymers are modeled by the Warner law, whereas 
biological molecules are described by the Marko-Siggia law. 
Relaxation times were computed by means of the variation- 
iteration method ^2^. For Cs = Cc they can be obtained by 
solving an eigenvalue problem for a Heun equation |15|]. 
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Figure 4: Bmnk-Koch-Lion flow: t-cd/r vs. Wi — Xt from 
Brownian Dynamics simulations; (a) PEO; (b) PAM. 

short-correlated flow. It is worth noting that the above 
definition provides an operational method to measure trd 
that can be implemented in experiments. 

In summary, we have shown that the equilibrium con- 
figuration of a polymer in a fiow, as well as the time a de- 
formed molecule takes on average to recover that config- 
uration, depend sensitively on the properties of the flow. 
In the vicinity of the coil-stretch transition the character- 
istic relaxation time is much longer than the Zimm time, 
both in elongational and random flows. In other words, 
the effective Wi differs considerably from the "bare" 
one [s^ . This effect is strongly amplified when the drag 
coefficient depends on the conformation of the polymer, 
and may play an important role in drag-reducing turbu- 
lent flows, where the strain rate often fluctuates around 
values typical of the coil-stretch transition "s^ . Our con- 
clusions thus suggest that the conformation-dependent 
drag should be included as a basic ingredient of contin- 
uum models of polymer solutions, calling for further the- 
oretical, experimental and numerical study. 
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